Specific heat of the S=| Heisenberg model on the kagome lattice: high-temperature 

series expansion analysis 
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We compute specific heat of the antiferromagnetic spin-i Heisenberg model on the kagome lattice. 
We use a recently introduced technique to analyze high-temperature series expansion based on the 
knowledge of high-temperature series expansions, the total entropy of the system and the low- 
temperature expected behavior of the specific heat as well as the ground-state energy. In the case 
of kagome-lattice antiferromagnet, this method predicts a low-temperature peak at T/J < 0.1. 
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PACS numbers: 75.10.Jm, 75.40.Cx 65.40.-b 65.40.Ba 

I. INTRODUCTION 

We consider the nearest-neighbor Heisenberg model on 
the kagome lattice: 



T~L — 2 S ' Si ■ Sj 



(1) 



Because of its unconventional 



properties, the spin-^ 
kagome antiferromagnet (KAF) has been subject to an 
intense activity these last years. All studies agree that 
this frustrated two-dimensional magnet 
ranged magnetic order at zero temperature.l 
Exact diagonalization studies have established that the 
low-energy spectrum of the kagome-lattice Heisenberg 
antiferromagnet has a large number oLspin-singlet states 
before the first spin 1 excited state.EJ Among the dif- 
ferent theories developed to explain this unconventional 
spectrum, short-range Resonating VaLgnce-Bond (RVB) 
pictures have been proposed.Mii3't3 : ESEj 

The high-temperature (HT) expansion of the specific 
heat as-been computed up to order 1/T 16 by Elstner and 
YoungJlj We have checked and extended this series to or- 
der 1/T 17 . The additional term for the specific heat per 
site is given by c v (T) = |/3 2 + - • •+ 18 Sg 3 /3 17 . Elst- 
ner and Young analyzed the series through conventional 
Pade approximants with the additional constraint that 
the specific heat must vanish at T = 0. At the highest 
orders, they found a specific heat curve with a single max- 
imum around T = 1.3 but with a large entropy deficit of 
about 40%: / °° c v {T)/TdT ~ 0.61n(2). They concluded 
the existence of a low-temperature structure correspond- 
ing to an entropy of about 40% of In 2 and claimed that 
this low-energy structure could not be accessed from the 
high-temperature expansion of the specific heat. They 
argued that even though the classical kagome antiferro- 
magnet has a non- vanishing ground-state entropy, quan- 
tum fluctuations in the spin-^ model are expected to lift 
this degeneracy. 



In this paper, we revisit the question of the specific 
heat with the help of a new method to analyze high- 
temperature series data. Comparevi|to the usual Pade 
approximant approach, this methodtZI takes advantage of 
additional information on the system: the two sum rules 
on the energy and on the entropy are exactly satisfied. 
In many simple systems (one- and two dimensional ferro- 
or antiferromagnets), this technique allows to compU|te| 
accurately the specific heat down to zero temperaturefH 
which is not the case if one does a direct Pade analysis of 
the series. For the present kagome model we show that 
this method provides rich semi-quantitative informations 
on the specific heat curve, although a full convergence 
down to zero temperature cannot be achieved. 



II. DIRECT HIGH-TEMPERATURE 
EXPANSION OF THE SPECIFIC HEAT 

We xeproduce here the first attempt by Elstner and 
YoungEEl to compute the specific heat from its high- 
temperature expansion alone. We use Pade approxi- 
mants to extrapolate the series. We impose the specific 
heat to vanish at low temperature as T, T 2 or T 3 . At 
orders 9 to 17, only 6 such approximants do not develop 
poles or zeros in the interval T g]0, oo] (Fig. Q). One 
should notice at this point that the remaining Pade ap- 
proximants agree reasonably well down to zero tempera- 
ture. This is usually not the case in 2D antiferromagnet 
where even the position of the peak (T ~ 1) can hardly be 
obtained by the use of direct Pade approximants to the 
series for the specific heat.E3 From this point of view, the 
HT series expansion of the kagome model seems to have 
a faster convergence than models such as the triangular- 
lattice antiferromagnet. 

By integration of these approximants, we evaluate the 
ground-state energy eo = J c v (T)dT and the ground- 
state entropy s = log(2) — J °° c v (T)dT. These values 
are indicated in Fig. The ground-state energy is about 
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FIG. 1: (color online) Specific heat obtained from Pade ap- 
proximants to the high-temperature series of c v (T) . Only ap- 
proximations of degrees [u, u+1], (resp. [it, u+2] and [u, it+3]) 
are considered. They vanish at zero temperature as T (resp. 
T 2 and T 3 ). The only 6 such approximants from order 9 to 
17 which are positive on the positive real axis are shown. The 
ground-state energy per site eo and entropy so obtained by 
integrating these Pade approximations are indicated. 



FIG. 2: (color online) Left: Series expansion for s(e) (entropy 
per site) as a function of the energy per site e. The result of 
the bare series are displayed for orders from /3 11 to f3 17 . The 
dashed vertical lines indicate upper (e = —0.84267) and lower 
(e = —0.909952) rigorous biWKids on the ground-state energy 
in the thermodynamic limit.t3 Right: Same data as on the left 
panel but plotted as a function of temperature T = l/s'(e). 
For each curve the lowest temperature corresponds to e = 
—0.909952 (lower bound). The black segment corresponds to 
e = -0.84267 (upper bound). 



—0.845, only slightly higher (0.02) than estimations ob- 
tained from exact diagonalizations. The entropy deficit- 
is very large: 0.3 (40% of log(2)). Elstncr and YoungH 
argued that a low-temperature peak should be present 
in the specific heat in order to compensate the deficit 
of 40% of log(2)£j However this low-temperature peak 
should "contain" almost no energy (2%), which means 
that such peak would have to occur at very low temper- 
atures. In order to estimate this temperature, one can 
add a (5-function peak to the curves of Fig. [l] in order to 
recover the correct energy and entropy variations. This 
constrains both the location (temperature T{) and the 
weight of the £-peak. By averaging over the different 
curves of Fig. one finds T x ^ 0.05 (resp. T x ^ 0.08) for 
a ground-state energy eo = —0.865 (resp. eo = —0.875). 
These estimates are in agreement with the conclusions of 
the more elaborate treatment described below. 



III. ENTROPY METHOD 

In this section we briefly summarize the method we use 
to compute the specific heat. More details can be found 
in Ref. The specific heat c v and the temperature T 
can be obtained from the entropy s as a function of the 
energy e using basic thermodynamic relations: 



T(e) = l/s'(e) 
s'(e) 2 



c„(e) 



s"(e) 



(2) 
(3) 



From Eq. M one can conver a high-temperature se- 
ries for c v (T — > oo) into a series for s(e — » 0) (e = at 
T = oo for the Hamiltonian of Eq. pT). The truncated 
series are plotted in Fig. |2| Using Eq7|^, the entropy can 



be plotted as a function of temperature (right of Fig. ||). 
A good convergence is observed down to relatively low 
energies (e ~ —0.75) but the corresponding entropy re- 
mains very large (more than 60% of ln(2)), although the 
ground-state energy is not much lower (the ground-state 
energy lies between the dashed vertical lines in Fig. ||). 
These result are consistent with a direct analysis of the 
series for c v {T) (Fig. |l|). In addition, it appears that the 
"true" s(e) must be bent downward below the curves of 
the truncated series (shown in Fig. |^) between eo and 
~ —0.75 in order to reach s = at e = e . Due to Eq. 0, 
this almost certainly implies a low-energy (and therefore 
low-temperature) peak in c v (T). This paper makes this 
idea more precise by computing the specific heat obtained 
by forcing the entropy to vanish at e = eo. 

The advantage of working on s(e) rather than c v (T) 
is that a two-point Pade interpolation can be used to 
set the ground-state energy and the total entropy of 
the system. However, s(e) is singular at e = eo (since 
T = l/s'(e) — * when e — > eo). For this reason one 
cannot directly approximate s(e) by a rational fraction 
(Pade approximant). If we assume that the specific-heat 
behaves as 



c v ~ (T/co) a 



(4) 



at low temperature (co has the dimension of an energy) 
and s(eo) = 0, s(e) behaves as: 

(o + 1)=tt /e-eo \ 
s(e -> e ) ~ ; ( — — ) (5) 



co 



The quantity 



G(e) 



s(e) 



1+1 /a 



e - e 



(6) 
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is then non-singular at e = eo and can be approximated 
by a Pade formBs The series for s(e — ► 0) must there- 
fore be transformed into a series for G(e — > 0) before 
Pade approximants can then be computed in the usual 
way. In what follows all the Pade approximants will 
be approximations to this function G(e). If no finite- 
temperature phase transition is expected, all approxi- 
mants where G(e) has a pole or a zero, or where s'(e) 
or s"(e) vanishes somewhere in the interval ]eo,0[ must 
be discarded. The remaining ones are called "physical" 
for brevity. 



A. Low-temperature behavior of c„(T) 

Unlike some simpler magnets where the nature p£ the 
ground-state and elementary excitations is known,EZI the 
qualitative behavior of the specific heat when T — > 
is unknown, although a ~ T 2 scenario has been 
proposed.Blia However, one of the striking facts about 
the model is the unusually high-density of states imme- 
diately above the ground-state.£3 From this it is natu- 
ral to expect gapless elementary excitations. If we as- 
sume quasi-particles with a dispersion relation ~ fc 7 
we get a specific heat c v ~ T a with a = D/j in space 
dimension D. The (many body) density of states is 
p(E + W) ~ exp [N(W/N) a ^ a+1 ^] where N is the sys- 
tem size (consequence of Eq. |^ with e — eo = W/N) . For 
an energy W of order one above the ground-state, a den- 
sity of states p ~ was observed in exact spectra 
up to N = 36 sitesJij If this indeed holds up to the ther- 
modynamic limit, it would imply a — (7 = 00) and an 
extensive entropy atjero temperature. This is unlikely 
in the present mode£3 but this result points to a rather 
flat dispersion relation of the excitations, probably with 
7 > 1. In the following we will .consider the two cases 
7 = 1 (a = 2) and 7 = 2 (a = 1).E3 



B. Ground-state energy and convergence of the 
different Pade approximants 

In principle, the method above requires the knowledge 
of the ground-state energy eo- If the value of eo is ex- 
act we expect the procedure to converge to the exact 
c v if the number of known terms in the HT series in- 
creases to infinity. This is in agreement with our experi- 
ence on solvable models (such as the spin-i XY chain for 
instancet3) where the full series as well as eo are known 
exactly. Inversely, wrong values of eo cannot lead to any 
convergence as the limiting c v would have to satisfy the 
HT series at all orders but would have a different energy 
sum rule. As a consequence, when eo differs from the 
true ground-state energy, the physical approximants gets 
fewer (and/or more scattered) when the order of expan- 
sion gets larger. Of course, the smaller the error on eo 
the longer series is needed to observe this departure from 
convergence. From this we assume that the existence of 



a larger number of physical approximant is an indication 
that eo (and a) is closer to the exact value. However, be- 
cause a limited number of terms of the series are known, 
this only provides a qualitative information and does not 
allow determine the energy completely. 

From exact diagonalizations on systems with, up to 
36 sites, eo was evaluated by Waldtmann et alEl to be 
e = -0.865 ± 0.015 (see also Refs. @,||,||). Variational 
calculations as well as rigorous bounds on eo will be dis- 
cussed in a separate paper .EjO 

The specific heat curve can be rather sensitive to the 
choice of eo- Since eo is not exactly known, it is important 
to perform scans in order to see how the specific-heat 
curve depends on eo- We observe that, for some choice of 
eo many Pade approximants at a given order give almost 
the same specific heat curve whereas some other choice of 
eo leads to some significant scattering in the specific heat 
curves. This can conveniently be seen, for instance, by 
looking at the value of the different Pade approximants 
at e = eo. Since G(eo) and Co (defined by Eq. |^) are 
simply related by 



G(e ) = 



a + 1 
1+1/c 



cq a 



(7) 



we plot Co (which has a direct physical meaning in term 
of c v (T — > 0)) in Fig. |as a function of eo for all phys- 
ical Pade approximants at order (3 16 and /3 17 (both for 
c v (T) ~ T and ~ T 2 at low temperature). It turns out 
that Co is representative of the full specific-heat curve in 
the sense that if two Pade approximants give "close" val- 
ues of Co (say a relative difference less than 10 -3 ), their 
corresponding specific-heat curves are similar (typical rel- 
ative difference of 10 -2 ) for all temperatures. This low- 
temperature coefficient Co is therefore a useful quantity 
to monitor how the c v (T) result depends on the choice 
of the degree of the Pade approximant J 3 !] In all cases the. 
"optimal" energy region is around eo — —0.88 ± 0.02.E1I 
We also observe a gradual shift of the optimal region to 
higher energies as the order of the series is increased. 
We analyzed this effect and performed several extrapo- 
lations to the infinite-order limit (data not shown). It is 
not clear, however, that this indirect method to deter- 
mine the ground-state energy is more accurate than the 
other available estimatesJla 



C. Low temperature peak in c v (T) 

The curves corresponding to all physical approximants 
at order (3 17 for e = —0.865, e = —0.88 and eo = 
—0.89, and for c v ~ T and c v ~ T 2 are shown in 
Fig. ^. Although some uncertainties remain concerning 
the ground-state energy of the model as well as the low- 
temperature behavior of the specific heat, the results are 
relatively well converged down to T ~ 0.7 and the loca- 
tion of the high temperature peak is almost independent 
from the unknowns, (fjp .and a) and is in agreement with 
previous studies BH'E3E3 In addition, all the scenarios we 
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FIG. 3: (color online) Zero temperature limit of Co (see 
Eq. ^|) for the different Pade approximants as a function of the 
ground-state energy. Top: ac„~ (T/co) behavior is assumed 
at low temperature. Bottom: c v ~ (T/co) 2 . The degree u of 
the numerator of each approximant is given. The degree of 
the denominator is d — n — u where n is the order of the series. 
Left: Order n — 16. Right: Order n — 17. The number of 
physical approximant (iVp a d c ) is plotted as a function of eo in 
each lower panel. 



investigated give rise to a low-temperature peak (or a 
shoulder) in the specific heat at T ~ 0.02 ~ 0.1. 

We also looked at the order dependence of the spe- 
cific heat curves. For a given value of the ground-state 
energy some approximants give similar curves for c v (T) 
while some other are "isolated". The later ones can be 
recognized-as isolated curves in Fig. ||. According to our 
experiencalJ with this method, those isolated approxi- 
mants do not reflect the convergence to the true func- 
tion. We obtained the results of Fig. || by keeping only 
the approximants which value of cq is at less than 3.10 -3 
from the c' of another approximant. This selection was 
repeated from orders /3 13 to /3 for the six combinations 
of ground-state energies and low-temperature behaviors 
used before. As one can see, the low-temperature struc- 
ture appears to be a robust feature, although a conver- 
gence of the full curve is not reached for T < 0.6. Still, 
a better convergence as a function of the order of the 
series (and a larger number of physical Pade approxi- 
mants) is obtained when the ground-state energy is low 
(eg = —0.89 or eq = —0.88). This suggests that the 




FIG. 4: (color online) Specific heat computed at order (3 17 
with the ground-state energies eg — —0.89 (top), eo = —0.88 
(center) and and eo = —0.865 (bottom). Left: c v ~ T. Right: 
Go ~ T 2 . The different curves correspond to all the physi- 
cal Pades approximants. The degree u of the numerator is 
indicated (the denominator has degree d = 17 — u). 



actual value of eo may be lower than —0.865, although 
eo = —0.89 is jjpabsdaly too low (compared to the avail- 
able estimatesBBElefl). 

For iV = 18 spins, exact diagonalizationsEl gave a low- 
temperature peak of the specific heat at T ~ 0.2 and 
c v ~ 0.17. An hybrid methodic based on exact diago- 
nalizations and high-temperature series expansion gave 
a peak at T ~ 0.2 and c v ~ 0.17 for N — 36 (see also 
Ref. @). Quantum Monte-Carlo simulations for N = 72 
spinsH indicated that a peak may exist below T ~ 0.3 for 
this system. Those results obtained for small systems are 
qualitatively consistent with those of Fig. ^| but our peak 
is located at a lower temperature by at least a factor of 
two. We think that this discrepancy is likely to be due 
to finite-size effects in previous studies. 
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IV. CONCLUSIONS 



By means of a detailed high-temperature series anal- 
ysis we provided quantitative estimates for the spe- 
cific heat curve of the spin-i Heisenberg antiferromag- 
net on the kagome lattice. Those results show a low- 
temperature peak in the specific heat of the model for 
T < 0.1, although its precise location cannot be deter- 
mined due to uncertainties on the ground-state energy. 
The corresponding degrees of freedom are also respon- 
sible for the large density of singlet states observed in 
exact diagonalization studies but their nature, as well 
as the nature of the ground state itself, remains to be 
explained. 
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FIG. 5: (color online) Specific heat curves from order /3 13 
to f3 17 with e = -0.89 (top), e = -0.88 (center) and 
eo = —0.865 (bottom). The degree of the numerator of 
each approximant is indicated and each column corresponds 
to a given order (/3 n ) of the series (n = 13, 14, 15, 16 and 
17 from left to right). Left panels: c v ~ T. Right panels: 
Cv ~ T 2 . In all cases the specific heat shows a maximum 
around T ~ 1.3 (corresponding to e(T) ~ —0.7, see Fig. ^) 
and a low-temperature peak (or shoulder). 
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J. J. FarnellpR- F. Bishop, 
Bernhard et alB used a mean-field decoupling on correla- 
tion functions to compute the specific heat of this system. 
Their result are qualitatively similar and also shows a large 
entropy deficit (46% of log(2)). 

Eq. ij is slightly different from the original choice made in 
Ref.]17|. 

See Ref. |l^ for a related model where this property is ex- 
plicitly demonstrated. This, however, requires some fine 
tuning of the Hamiltonian and there is no reason to think 
that the kagome Heisenberg model may have such very 
special property. From finite-size spectra it can however 
be difficult to resolve the difference between a model with 
an extensive entropy at T = and a system with vanishing 
entropy at T = but a large low-temperature peak in the 
specific heat for T <C J. We stress that the T = entropy 
that follows from p ~ 1.15 is consistent with that of the 
low-temperature peak predicted by the present method: 
about 20%-25% of log(2). 

The possibility of having a gapped spectrum with a (very) 
small gap cannot be excluded. This would be the case if 
the system realizes a Z2 spin-liquid, as predicted in large- 
N approaches^ or if thpnW-o«*id-state has some valence- 
bond long-ranged order .EjO'EiJ A thermally activated be- 
havior Cv ~ f»sp(— A/T) can be treated within the en- 
tropy methodtj and we applied it to the present model. 
The gaps we obtained are rather large (of order 1) but the 



specific-heat curves also exhibited a low-temperature peak 
consistent with the c„ ~ T a results discussed below. 
The expectation value of the Hamiltonian in any first- 
neighbor valence-bond covering is e = —0.75. By opti- 
mizing the wave-function in the vicinity of each triangle 
without any dimer, Elser obtained a variational state with 
energy e = — 0t8333. This rigorous upper bound was re- 
fined by PierreE3 who used an improved variational state 
to prove that eo < —0.84267. And exact lower bound 
can be obtained by remarking that the ground-state en- 
ergy per triangle in the kagome AF is necessarily larger 
than the ground-state energy of an isolated triangle. This 
gives eo > —1. By extending this reasoning to larger clus- 
ters we can improve this lower bound. We can prove that 
eo > —0.909952. by exactly diagonalizing a 24-site clus- 
ter with open boundary conditions. Farnell, et al. used 
a coupled cluster method and predicted en — — 0.8504.E3 
They recently improved their calculations and obtained 
e ~ -0.86208 by fitting their n th -order results by 1/n 2 
corrections. We observed that their bare data are bet- 
ter described by 1/n corrections. In that case a fit gives 
e ~ -0.875. 

The physical meaning of this observation is that if two 
specific heat curves have i) the high-temperature expansion 
up to some relatively high order ii) the same ground-state 
energy, iii) the same entropy and iv) the same c v ~ (T/co) a 
limit at low temperature they must be very "similar" . 
For lower energies the curves gets closer but the number 
of physical Pade approximants actually gets smaller and 
smaller, as can be seen in the lower panels of Fig. pi 



